home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2001 May / SGI IRIX Base Documentation 2001 May.iso / usr / share / catman / p_man / cat3 / complib / CINVIT.z / CINVIT
Encoding:
Text File  |  1998-10-30  |  6.1 KB  |  133 lines
  1.  
  2.  
  3.  
  4. ____CCCCIIIINNNNVVVVIIIITTTT((((3333FFFF))))                                                        ____CCCCIIIINNNNVVVVIIIITTTT((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CINVIT, SCINVIT  -  EISPACK routine.  This subroutine finds those
  10.      eigenvectors of A COMPLEX UPPER Hessenberg matrix corresponding to
  11.      specified eigenvalues, using inverse iteration.
  12.  
  13.  
  14. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  15.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee  cccciiiinnnnvvvviiiitttt((((nnnnmmmm,,,, nnnn,,,, aaaarrrr,,,, aaaaiiii,,,, wwwwrrrr,,,, wwwwiiii,,,, sssseeeelllleeeecccctttt,,,, mmmmmmmm,,,, mmmm,,,, zzzzrrrr,,,, zzzziiii,,,,
  16.           1            ierr, rm1, rm2, rv1, rv2)
  17.           iiiinnnntttteeeeggggeeeerrrr          nnnnmmmm,,,, nnnn,,,, mmmmmmmm,,,, mmmm,,,, iiiieeeerrrrrrrr
  18.           ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn aaaarrrr((((nnnnmmmm,,,,nnnn)))),,,,  aaaaiiii((((nnnnmmmm,,,,nnnn)))),,,, wwwwrrrr((((nnnn)))),,,, wwwwiiii((((nnnn))))
  19.           ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn zzzzrrrr((((nnnnmmmm,,,,mmmmmmmm)))),,,, zzzziiii((((nnnnmmmm,,,,mmmmmmmm))))
  20.           ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn rrrrmmmm1111((((nnnn,,,,nnnn)))),,,,  rrrrmmmm2222((((nnnn,,,,nnnn)))),,,, rrrrvvvv1111((((nnnn)))),,,, rrrrvvvv2222((((nnnn))))
  21.           llllooooggggiiiiccccaaaallll          sssseeeelllleeeecccctttt((((nnnn))))
  22.  
  23.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee sssscccciiiinnnnvvvviiiitttt((((nnnnmmmm,,,, nnnn,,,, aaaarrrr,,,, aaaaiiii,,,, wwwwrrrr,,,, wwwwiiii,,,, sssseeeelllleeeecccctttt,,,, mmmmmmmm,,,, mmmm,,,, zzzzrrrr,,,, zzzziiii,,,,
  24.           1            ierr, rm1, rm2, rv1, rv2)
  25.           iiiinnnntttteeeeggggeeeerrrr          nnnnmmmm,,,, nnnn,,,, mmmmmmmm,,,, mmmm,,,, iiiieeeerrrrrrrr
  26.           rrrreeeeaaaallll             aaaarrrr((((nnnnmmmm,,,,nnnn)))),,,,  aaaaiiii((((nnnnmmmm,,,,nnnn)))),,,, wwwwrrrr((((nnnn)))),,,, wwwwiiii((((nnnn))))
  27.           rrrreeeeaaaallll             zzzzrrrr((((nnnnmmmm,,,,mmmmmmmm)))),,,, zzzziiii((((nnnnmmmm,,,,mmmmmmmm))))
  28.           rrrreeeeaaaallll             rrrrmmmm1111((((nnnn,,,,nnnn)))),,,,  rrrrmmmm2222((((nnnn,,,,nnnn)))),,,, rrrrvvvv1111((((nnnn)))),,,, rrrrvvvv2222((((nnnn))))
  29.           llllooooggggiiiiccccaaaallll          sssseeeelllleeeecccctttt((((nnnn))))
  30.  
  31.  
  32.  
  33. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  34.      On INPUT
  35.  
  36.      NNNNMMMM must be set to the row dimension of two-dimensional array parameters
  37.      as declared in the calling program dimension statement.
  38.  
  39.      NNNN is the order of the matrix.
  40.  
  41.      AAAARRRR and AI contain the real and imaginary parts, respectively, of the
  42.      Hessenberg matrix.
  43.  
  44.      WWWWRRRR and WI contain the real and imaginary parts, respectively, of the
  45.      eigenvalues of the matrix.  The eigenvalues must be stored in a manner
  46.      identical to that of subroutine  COMLR, which recognizes possible
  47.      splitting of the matrix.
  48.  
  49.      SSSSEEEELLLLEEEECCCCTTTT specifies the eigenvectors to be found.  The eigenvector
  50.      corresponding to the J-th eigenvalue is specified by setting SELECT(J) to
  51.      .TRUE.
  52.  
  53.      MMMMMMMM should be set to an upper bound for the number of eigenvectors to be
  54.      found.  On OUTPUT
  55.  
  56.      AAAARRRR , AI, WI, and SELECT are unaltered.
  57.  
  58.      WWWWRRRR may have been altered since close eigenvalues are perturbed slightly
  59.      in searching for independent eigenvectors.
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. ____CCCCIIIINNNNVVVVIIIITTTT((((3333FFFF))))                                                        ____CCCCIIIINNNNVVVVIIIITTTT((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      MMMM is the number of eigenvectors actually found.
  75.  
  76.      ZZZZRRRR and ZI contain the real and imaginary parts, respectively, of the
  77.      eigenvectors.  The eigenvectors are normalized so that the component of
  78.      largest magnitude is 1.  Any vector which fails the acceptance test is
  79.      set to zero.
  80.  
  81.      IIIIEEEERRRRRRRR is set to Zero       for normal return, -(2*N+1)   if more than MM
  82.      eigenvectors have been specified, -K         if the iteration
  83.      corresponding to the K-th
  84.         value fails, -(N+K)     if both error situations occur.
  85.  
  86.      RRRRMMMM1111 , RM2, RV1, and RV2 are temporary storage arrays. The ALGOL procedure
  87.      GUESSVEC appears in CINVIT in line.  Calls PYTHAG(A,B) for sqrt(A**2 +
  88.      B**2).  Calls CDIV for complex division.  Questions and comments should
  89.      be directed to B. S. Garbow, APPLIED MATHEMATICS DIVISION, ARGONNE
  90.      NATIONAL LABORATORY
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.